Structural and dielectric properties of amorphous Zr02 and Hf02 
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Zirconia (Zr02) and hafnia (Hf02) are leading candidates for replacing Si02 as the gate insulator 
in CMOS technology. Amorphous versions of these materials (a-Zr02 and a-Hf02) can be grown 
as metastable phases on top of a sihcon buffer; while they tend to recrystallize during subsequent 
annealing steps, they would otherwise be of considerable interest because of the promise they hold for 
improved uniformity and electrical passivity. In this work, we report our theoretical studies of a-Zr02 
and a-Hf02 by first-principles density-functional methods. We construct realistic amorphous models 
using the "activation-relaxation" technique (ART) of Barkema and Mousseau. The structural, 
vibrational, and dielectric properties of the resulting models are analyzed in detail. The overall 
average dielectric constant is computed and found to be comparable to that of the monoclinic 
phase. 
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PACS numbers: 77.22.-d, 61.43.Bn, 63.50.+X, 71.23.Cq 
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I. INTRODUCTION 

In recent years, a major thrust of applied development 
in the semiconductor industry has been the search for 
materials that could replace Si02 as the gate dielectric 
in CMOS technology. Conventional scaling of Si02 would 
require a gate-dielectric thickness shrinking below 1 nm 
in the near future, according to the International Tech- 
nology Road Map for Semiconductors.^ In this regime, 
conventional thermally-grown Si02 is expected to fail be- 
cause of issues of tunneling leakage current and reliability. 
There are therefore strong incentives to identify materials 
with dielectric constant e (or K) much larger than that 
of Si02 (at e=3.9), as these could be grown as thicker 
films while still providing the needed capacitance. 

Among the most promising materials are Zr02 and 
Hf02, which do have much higher dielectric constants 
and some other positive features as well (e.g., chemical 
stability). However, one of the great advantages of Si02 
has been the fact it forms an amorphous oxide (a-Si02), 
thus allowing it to conform to the substrate with enough 
freedom to eliminate most electrical defects at the inter- 
face. On the other hand, Zr02 and IIf02 (which are very 
similar in many of their physical and chemical proper- 
ties) are refractory materials, their melting temperature 
being 2988 K and 3085 K respectively. They are not good 
glass- formers; a-Zr02, for example, has been shown to 
recrystallize during growth if the growth temperature is 
too high (with tetragonal and monoclinic phases starting 
to appear above ^500°C and '^700°C respectively) pi 
Thus, while Zr02 and Hf02 can be grown as metastable 
amorphous phases on Si using low-temperature deposi- 
tion techniques, films of this type unfortunately tend to 
recrystallize during the subsequent annealing steps that 
are required in current industrial fabrication processes. 

Nevertheless, it is possible that admixing (alloying) 
with Si, Al, N, or other chemical constituents, or other 
strategies yet to be identified, may help to mitigate the 
recrystallization problem and stabilize the amorphous 



phasei^ In any case, it may be advantageous to study 
amorphous structures as a first step in understanding 
why the recrystallization is facilitated. With these moti- 
vations, we have embarked on a theoretical study of the 
structural and dielectric properties of amorphous zirco- 
nia and hafnia (a-Zr02 and a-Hf02). In earlier work^ 
we generated models of a-Zr02 using an ab-initio molec- 
ular dynamics approach in a plane-wave pseudopotential 
framework, and studied their structural and dielectric 
properties. In a subsequent publication,^ we described an 
alternate and more efficient method of generating amor- 
phous structures based on the activation-relaxation tech- 
nique (ART) of Barkema and Mousseau,^ and presented 
a few preliminary results on amorphous Zr02 obtained 
using this approach. These calculations were carried 
out using a local-orbital basis approach embodied in the 
SIESTA code packagei^ In the present work, we system- 
atically use SIESTA for all of the reported calculations. 
We first briefly review our implementation of the ART 
approach, and then embark on a systematic description 
of the generated a-Zr02 and a-IIf02 structures and a 
comparison of their computed electronic, structural, and 
dielectric properties. 

The paper is organized as follows. In Sees. II and III 
we describe the details of our electronic-structure calcu- 
lations and of our procedure for generating amorphous 
models, respectively. We then describe the results of our 
calculations for Zr02 in Sec. IV and for Hf02 in Sec. V. 
We also present some results concerning the Born dynam- 
ical effective charges and the dielectric activity activity 
for both systems in Sec. VI. We then finish with a con- 
clusion in Sec. VII. 



II. COMPUTATIONAL DETAILS 

We performed ab-initio density-functional theory 
(DFT) calculations! in the local-density approximation 
(LDA)^ using a local-orbital expansion of the Bloch wave- 
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function as implemented in the SIESTA code.- The use of 
locahzcd-orbital basis allows the retention of good accu- 
racy in the DFT calculations at a reduced computational 
cost relative to plane-wave codes. A cutoff of 150 Ry 
was used for the expansion of the charge density, and all 
calculations were performed using F-point Brillouin zone 
sampling only. The dynamical matrix and Born effective 
charges were computed using the SIESTA via a finite- 
difference approach in which each atom was displaced by 
±0.05 A in the three cartesian directions and the forces 
and polarization changes were computed. 

Two kind of basis sets were employed in the calcu- 
lations. A minimal single-C basis set was used for the 
amorphization procedure, which required a large number 
of structural relaxations, while a larger triple-^ basis set 
was used for the structural relaxation of the final struc- 
ture and for the calculation of the lattice dielectric con- 
stant and of the infrared activity. Such a two-stage pro- 
cedure (amorphization with a lower-level theory followed 
by annealing and relaxation with a higher-level theory) is 
a natural way of establishing a trade off between compu- 
tational cost and accuracy; it has many precedents in the 
literature, as for example the generation of an amorphous 
structure by the use of empirical potentials followed by a 
first-principle annealing run in the work of Sarnthein et 

aa 

Because we were concerned that the minimal basis 
might not accurately reproduce the delicate interplay be- 
tween ionic and covalent bonding in Zr02 and IIf02, the 
parameters for the minimal basis set were optimized fol- 
lowing the prescription of Ref. |lo| in order to ensure the 
correct energy ordering of the cubic, orthorhombic and 
monoclinic phases of Zr02 and Hf02. We then obtain 
Ecnb - Etct = 62 meV and Etet - E^ono = 89 meV for 
Zr02, to be compared with values of 44 meV and 45 meV, 
respectively, computed using plane- wave methods and re- 
ported in Ref. Here £^cub, -E^tet, and i^mono are the 
energies per formula unit of the cubic, tetragonal, and 
monoclinic phases, respectively of Zr02. Similar results 
are found for IIf02. The relaxed structural parameters 
are also found to agree well. Thus, despite their sim- 
plicity, our optimized minimal basis sets for Zr02 and 
Hf02 are able to capture the essentials of the structural 
energetics of the three crystalline phases. 



III. AMORPHIZATION PROCEDURE 

In a previous studyj^ we carried out ab-initio molecular 
dynamics (MD) simulations of a 96-atom supercell in a 
melt-and-quench fashion in order to generate a structural 
model for a-Zr02. Because of the short time interval ac- 
cessible to simulations, the cooling rate (3.4x10^'* K/s) 
was far beyond the fastest cooling rate that can be ob- 
tained experimentally by pulsed laser techniques. The 
fast cooling rate does not allow long-time-scale relaxation 
of system, which might be important in the case of rela- 
tively poor glass formers such as Zr02 and IIf02. 



In order to see if the fast cooling rate could bias the 
resulting amorphous structure, we generated an indepen- 
dent amorphous sample^ by using the ART event-based 
structural-evolution approach.^ An event-based method 
leads to an accelerated dynamics and to a better sam- 
pling of long-time-scale modes in glasses. Rather than 
following the irrelevant details of the atomic motions as 
atoms vibrate back and forth about their average po- 
sitions for long periods between activated events, ART 
focuses on simulating jumps over the barriers that sep- 
arate the different basins of attraction of different local 
minima in the energy landscape. Thus, it requires a given 
computational effort per activated hop instead of per vi- 
brational period. For disordered systems, which tend to 
show slow evolution, these two time scales may differ by 
many orders of magnitude, and this event-based tech- 
nique therefore allows a much faster simulation. 

We implemented the ART method as a driver for the 
ab-initio code SIESTA J- The flowchart of our implemen- 
tation of the ART method can be found in Ref. 0. As 
reported there, we found that, starting from a 96-atom 
supercell made from the perfect cubic fluorite crystal 
structure, a number of Monte-Carlo (MC) trials equal to 
5 times the number of atoms (i.e., ^500) was sufficient 
to produce a good amorphous structure. In our case, a 
MC temperature of 3000 K produced an acceptance ra- 
tio of 13%. Moreover, we found that during the first 
50 MC trials, the acceptance ratio was higher (~ 40%) 
and dropped down in the subsequent trials. We did not 
attempt to perform longer simulations at a lower MC 
temperature in order to "anneal" the system further. 



IV. RESULTS FOR ZrOa 

An important issue is the density of the amorphous 
phase, which, to our knowledge, is not accurately known 
experimentally. This issue was explored in our previous 
papers^ and we found that among the initial guesses, a 
density between 4.86 g/cm^ and 5.32 g/cm^ generates a 
robust amorphous structure for Zr02;^^ 

The ART simulation was performed at a constant vol- 
ume corresponding to a density of 5.32 g/cm'^. A snap- 
shot of the system is essentially indistinguishable, at 
the visual level, from that of the melt-and-quench MD- 
generated system of Ref. 3. The corresponding distribu- 
tion of coordination numbers is shown in Fig.^ We find a 
prevalence of 7- and 6-coordinated Zr over 8-coordinated 
Zr and a prevalence of 3- and 4-coordinated oxygens, sim- 
ilar to what was found for the MD-generated structure. 
(For comparison, recall that the monoclinic structure has 
7- fold cations and equal numbers of 3- and 4- fold anions.) 

One may wonder if such a spread of coordination num- 
bers is actually due to intrinsic coordination defects such 
as dangling bonds or non-bridging oxygens (which are of- 
ten present in other oxides such as Si02). We monitored 
the electronic density of states at every MC accepted 
move and we found that occasionally a defective struc- 
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FIG. 1: a-Zr02 distribution of coordination numbers resulting 
from the activation-relaxation (ART) simulation. Zr and O 
atoms are indicated by filled and open bars, respectively. 



FIG. 3; a-Hf02 distribution of coordination numbers result- 
ing from the activation-relaxation (ART) simulation. Hf and 
O atoms are indicated by filled and open bars, respectively. 



ture (i.e., with some levels in the gap) is accepted. We 
also found, however, that the coordination defects are 
promptly saturated in subsequent MC moves. The final 
amorphous sample is a good insulator with an electronic 
gap of ~ 3.4 eV. 

Furthermore, the calculated phonon spectrum, shown 
in Fig. 12 is found to extend over approximately the same 
range of frequencies (50-800cm~^) and to show features 
similar to those of the sample obtained by the melt-and- 
quench MD simulation. The computed Born effective 
charges are slightly larger on average, Z*(Zr)=4-5.08 and 
Z*(0)=— 2.54. Our lattice dielectric tensor is 



17.9 -0.7 0.2 
eiatt = I -0.7 17.5 -0.5 
0.2 -0.5 14.1 



(1) 



yielding an average dielectric constant of 16.5, compared 
to the value of 17.6 of the MD-generated model, due 
to the fact that the vibrational spectrum is shifted to 
slightly higher frequencies. Assuming a value of too = 4.6 



for the high-frequency dielectric constant)^ the static di- 
electric constant yields a value of ~21, in good agree- 
ment with experimental valuesi^ and with the previous 
calculationi^ 



V. RESULTS FOR HfOa 

The amorphization of Hf02 was carried out in a man- 
ner very similar to that used for its twin oxide Zr02. 
For Hf02, however, we avoided the tedious procedure of 
performing independent melt- and- quench simulations in 
order to find a suitable density that can sustain an amor- 
phous structure. Instead, we started the ART simulation 
at the same volume as for Zr02, 38.45 A'^ per unit formula 
(corresponding to an Hf02 mass density of 9.09g/cm'^). 
Then, every 10 accepted MC moves, we performed a re- 
laxation of the supercell volume. The final amorphous 
Hf02 sample had a density of 9.39 g/cm'^ (corresponding 
to 37.22 A3 per unit formula, 3.2 % smaller than that of 
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FIG. 2: (a) a-Zr02 phonon density of states (DOS) vs fre- 
quency, (b) a-Zr02 infrared activity (phonon DOS weighted 
by Z\/ui1) vs frequency. 
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FIG. 4; (a) o-IIf02 phonon density of states (DOS) vs fre- 
quency, (b) a-Hf02 infrared activity (phonon DOS weighted 
by Z\/uj\) vs frequency. 
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FIG. 5: (a) Scatter plot of isotropically averaged atomic Z* 
values (vertical axis) vs atom type and coordination number 
(horizontal axis) for a-Zr02. Circles and diamonds denote 
O and Zr atoms, respectively, (b) Same but with "dielectric 
activity" (see Eq. (5) of Ref. |^ plotted vertically. 



FIG. 6: (a) Scatter plot of isotropically averaged atomic Z* 
values (vertical axis) vs atom type and coordination number 
(horizontal axis) for a-Hf02. Circles and diamonds denote 
O and Hf atoms, respectively, (b) Same but with "dielectric 
activity" (see Eq. (5) of Ref. 0) plotted vertically. 



Zr02). Experimentally, the unit cell volumes of the mon- 
oclinic and cubic phases are ~2% larger in Hf02 than in 
the corresponding Zr02 phasesii 

The distribution of coordination numbers of the re- 
sulting a-Hf02 sample is shown in Fig. |21 The aver- 
age coordination number shows a prevalence of 7- and 
6-coordinated Hf over 8-coordinated Hf, and a preva- 
lence of 3- and 4-coordinated oxygens, as in the mon- 
oclinic structure. Compared to a-ZiO^ (Fig.P) there is a 
slightly higher prevalence of 3-coordinated oxygens, while 
the distribution of 6- and 7-coordinates cations is broadly 
similar. However, it is doubtful whether the small differ- 
ences between Zr02 and Hf02 that are visible in Figs. 
and 131 are statistically significant. 

The generated a-Hf02 sample is found to be insulat- 
ing, with a gap of ~3.8eV. Its calculated phonon spec- 
trum, shown in Fig. 01 extends over approximately the 
same range of frequencies (50-800 cm^^) as for a-Zr02- 
Compared to Fig. [2Ib), the infrared activity for Hf02 in 
Fig-EIb) shows a broader feature in the 50-400 cm^^ fre- 
quency range. The computed Born effective charges are 
found to be smaller, on average, compared to those of 
a-Zr02: Z*(Hf)=-)-4.8 and Z*(0)=-2.4. Our computed 
lattice dielectric tensor 

/ 23.4 -5.1 -1.2 \ 
eiatt = -5.1 14.4 0.2 (2) 
y -1.2 0.2 12.6 ) 

yields an average lattice dielectric constant of 16.8. As- 
suming a value of ~ 5 for the high-frequency dielectric 
constant [eoo)r^ the static dielectric constant yields a 
value of ~ 22, confirming again the striking similarity be- 
tween a-Zr02 and a-Hf02. 



VI. BORN CHARGES AND DIELECTRIC 
ACTIVITY 

Clearly it is desirable to understand more fully the 
various contributions to the lattice dielectric response of 
the amorphous forms of Zr02 and Hf02 . To this end, we 
decomposed various lattice properties by "atom type" 
(that is, by chemical species and coordination number) 
in the hope that such an analysis may provide further 
insight into our numerical results. This type of analysis 
was explained in detail and applied to the a-Zr02 sample 
generated using the MD melt-and-quench approach in 
Ref.H 

Our results for the ART-generated samples are pre- 
sented in Figs. |31 and Fig. ^ In the top panel of each 
figure, we report the Born effective charges as a function 
of the chemical species and coordination number. In both 
cases, the Born effective charges tend to be surprisingly 
independent of coordination number, and quite similar 
to the values found in the crystalline phasesiii However, 
a mild tendency of the effective charge to increase with 
increasing coordination number of the cation is visible 
especially for Zr; it is also present for Hf but is partially 
obscured by larger fluctuations of the effective charges in 
that case. 

The atom-resolved dielectric activity, defined in Eq. (5) 
of Ref. 0, is presented in Figs-ISfb) and|Bfb). This is es- 
sentially a measure of the contribution of each type of 
atom to the lattice dielectric constant. Some trends are 
visible, such as a stronger contribution by 3-fold than 
by 4- fold oxygens, and a tendency to have larger contri- 
butions from higher-coordinated cations. However, even 
within any one atom type, there is a surprising degree 
of variation, with some atoms contributing strongly and 
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others contributing very little. Since the effective charges 
do not have nearly such a large variation, we infer that 
the differences must arise because of differences in the 
force-constant matrix elements. That is, some atoms 
may be regarded as being strongly anchored in place so 
that they contribute only weakly to the dielectric activ- 
ity, while others participate in one or more soft modes 
and contribute strongly. 

VII. CONCLUSIONS 

In conclusion, we have applied ab-initio electronic 
structure methods to study the lattice dynamics and 
dielectric properties of amorphous high-if materials a- 
Zr02 and a-Hf02. We used the ART event-based tech- 
nique to generate structural models of these amorphous 
materials. These 96-atom supercell models display a dis- 
tribution of coordination numbers (mostly 3- and 4-fold 
for oxygen and 6-, 7- and 8-fold for Zr and Hf) but 
nonetheless remain insulating with a robust gap. The 
full force-constant matrix was computed for each super- 
cell model, and the phonon density of states was calcu- 
lated. In addition, the full Born-charge tensors were also 
computed for each model, and when combined with the 
force-constant information, the lattice dielectric response 
was obtained. The total dielectric constant is computed 
to be ~ 22 for both materials, comparable to that of the 
monoclinic phase. Finally, the Born charges and the con- 
tributions to the dielectric activity were further analyzed 
by atom type. The effective charges are relatively uni- 
form and are roughly similar to those of the crystalline 



phases, whereas there are strong variations in the indi- 
vidual atomic contributions to the dielectric activity re- 
sulting from variations in soft-mode participation. 

Several caveats are in order. First, while our ART 
procedure can be regarded as corresponding to a slower 
quench than in the previous MD-based calculations^^ it 
is still very fast relative to any experimental quench- 
ing time scale, so that our models might have more de- 
fects than the experimental systems of interest. Sec- 
ond, real amorphous systems are generally grown at suf- 
ficiently low temperature to avoid crystallization, and 
not quenched from a melt. Thus, there could be signifi- 
cant differences in the resulting structures (for example, 
a greater tendency to void formation in the experimental 
samples). Third, our samples are perfectly stoichiomet- 
ric and impurity- free, which will not generally be true of 
experimental samples. Finally, while a statistical analy- 
sis of an ensemble of ART-generated samples would be 
very desirable, the computational burden needed to carry 
out the above analysis for even a single sample is quite 
demanding. Nevertheless, we believe that the present re- 
sults reveal the broad qualitative features to be expected 
for these amorphous phases, and constitute an impor- 
tant step forward on the road to a better understanding 
of these important materials. 
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